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Abstract. We apply the analytic-numerical method of Roberts to determine the linear 
stability of time-reversible periodic simultaneous binary collision orbits in the symmetric 
collinear four body problem with masses 1, m, m , 1, and also in a symmetric planar 
four-body problem with equal masses. For the collinear problem, this verifies the earlier 
numerical results of Sweatman for linear stability. 



1. Introduction 

Recently, Roberts [9| desribed an analytic-numerical method for determining the linear 
stability of a symmetric periodic orbit of a Hamiltonian system. He applied this method to 
the time-reversible collision-free figure-eight orbit in the equal mass three body problem 
numerically discovered by Moore IS) and whose existence was proven by Chenciner and 
Montgomery IJ]- (Other such choreographic solutions were found numerically by Simd 
ifTDl ). Roberts' method shows that the figure eight orbit is linearly stable. The method 
uses the symmetries to factor a matrix similar to the monodromy matrix for the periodic 
orbit into an integer power of the product of two involutions. One of the two involutions 
depends on the linearized dynamics along only a part of the periodic orbit. For the figure 
eight this part is one-tweltfth of the full orbit since it has a symmetry group isomorphic to 
the group D3 x Z2 of order 12. (Here the dihedral group Dk is the group of symmetries of 
the regular ^-gon.) The eigenvalues of the product of the two involutions are then reduced 
to the numerical computation of a few real numbers. 

Schubart IliO) numerically discovered a singular periodic orbit in the collinear equal 
mass three-body problem. The orbit alternates between binary collisions. Henon |[3l ex- 
tended Schubart's numerical investigations to the case of unequal masses. Only recently 
did Venturelli lfT4l and Moeckel 13 prove the existence of the Schubart orbit when the outer 
masses are equal and the inner mass is arbitrary. The linear stability of the Schubart orbit 
was determined numerically by Hietarinta and Mikkola |4| revealing that linear stability 
occurs for some but not all of the choices of the three masses. Sweatman ( lfT2l and ITSl ) 
numerically found and determined the linear stability of a Schubart-like orbit in the sym- 
metric collinear four body problem with masses 1 , m, m, and 1 . This Schubart-like periodic 
orbit alternates between simultaneous binary collisions (SBC) and inner binary collisions. 
Ouyang and Yan fj] proved the existence of this orbit. In the regularized setting, this pe- 
riodic orbit has a symmetry group isomorphic to D2, of which both of the generators are 
time-reversing symmetries. Ouyang, Yan, and Simmons ID numerically found and then 
proved the existence of a singular periodic orbit in a symmetric planar four-body problem 
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with equal masses in which the four bodies alternate between different simultaneous binary 
collisions. In the regularized setting, this periodic orbit has a symmetry group isomorphic 
to D4, of which one of the generators is a time-reversing symmetry. The regularization of 
these singular periodic orbits is achieved by a generalized Levi-Civita type transformation 
and an appropriate scaling of time, as adapted from Aarseth and Zare |[T|. 

In this paper we apply the method of Roberts to prove the linear stability of the Schubart- 
like orbit in the symmetric collinear four body \,m,m,\ problem for certain values of m, 
and of the singular periodic orbit in the symmetric planar equal mass problem. In both 
settings, the linear stability is determined for the regularized equations only and is reduced 
to the rigorous numerical computation of a single real number Our linear stability analysis 
determines values of m in the interval [0, 50] in the collinear problem for which the singular 
periodic orbit is linear stable, and also shows that the 2D singular periodic orbit is linear 
stable. These examples support and extend the conjecture made by Roberts 19| that the 
only linearly stable periodic orbits in the equal mass «-body problem are those that exhibit 
a time-reversing symmetry. 

Our linear stability analysis confirms Sweatman's linear stability analysis lfT3l for the 
singular periodic orbit in the collinear four-body problem. Sweatman used a numerical 
perturbation technique to assess the stability of the singular periodic orbit when the masses 
are arranged from left to right as m\, 1112, 1112, and mi with the condition that m\ +1112 = 
2. Our mass parameter m is related to his mass parameter mi by m = (2 — mi) /mi. In 
terms of our mass parameter m, Sweatman's numerical results indicate that linear stability 
occurs when the value of m is smaller than approximately 2.83 and when it is larger than 
approximately 35.4, and is linearly unstable otherwise. 



2. Linear Stability of Periodic Orbits 

For a smooth function F defined on an open subset of M?", suppose that y{s) is a T- 
periodic solution of a Hamiltonian system z' = JDr{z) where ' ^ d/ds, 



J 



/ 

-/ 



and / is the appropriately sized identity matrix. The fundamental matrix solution X{s) of 
the linearized equations along y(i), 

(1) ^'=jD'msm, m=i 

is symplectic and satisfies X{s + T) = X{s)X{T) for all s. The matrix X{T) is commonly 
called the monodromy matrix for y, and it measures the non-periodicity of solutions to the 
linearized equations. The eigenvalues of X{T) are the characteristic multipliers of y, and 
determine the linear stability of the periodic solution y. Linear stability therefore requires 
that all of the multipliers lie on the unit circle. 

The characteristic multipliers may be obtained by solving ([T]i with different initial con- 
ditions. For an invertible matrix Yo, let Y{s) be the fundamental matrix solution to 

(2) ^'^JD^T{y{sm, m^Yo. 

By definition ofX{s), we know that Y{s) = X{s)Yo, and so X{T) = Y{T)YqK It follows 
that the matrix Yq^Y{T) is similar to the monodromy matrix i.e., 

X{T)^Y{T)Y,-'^Yo{Y,-'Y{T))Y,-K 

Thus the eigenvalues of Yq^Y{T) are identical to the characteristic multipUers. 
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2. 1 . Stability reduction using symmetry. The monodromy matrix for a periodic solution 
with special types of symmetry can be factored using some linear algebra and standard 
techniques in differential equations. We begin by reviewing the relevant factorization and 
reduction theory that are applicable to a wide range of symmetric periodic orbits commonly 
found in Hamiltonian systems. Proofs of the following statements can be found in ||9l . 

Lemma 2.1. Suppose that j(s) is a symmetric T —periodic solution of a Hamiltonian sys- 
tem with Hamiltonian Y and symmetry matrix S such that: 

(1) for some positive integer N, y{s + T /N) = Sj{s) for all s; 

(2) r(5z)=r(z); 

(3) SJ = JS; 

(4) S is orthogonal. 

Then the fundamental matrix solution X{s) to the linearization problem in ([T]i satisfies 

X{s + T/N)=SX{s)S'^X{T/N). 

Here of course, the notation means the transpose of S. We mention this because we 
are using the letter T in two distinct ways. 

Corollary 2.2. Given the hypothesis of Lemma \2.1\ the fundamental matrix solution X{s) 
satisfies 

X{kT/N)^S\s'^X{T/N)Y 

for any k CzN. 

A remark here is that if Y{s) is the fundamental matrix solution to Equation (|2|, then 

for any e N, the matrix Y{kT /N) factors as 

Y{kT/N)=S%(YQ^S^Y{T/N))''. 

Lemma 2.3. Suppose that y(s) is a T— periodic solution of a Hamiltonian system with 
Hamiltonian Y and time- reversing symmetry S such that: 

(1) for some positive integerN, '•/{—s + T/N) — Sj{s) for all s; 

(2) r(5z)=r(z); 

(3) SJ = -JS; 

(4) S is orthogonal. 

Then the fundamental matrix solution X{s) to the linearization problem in ([U satisfies 

X{-s + T/N)=SX{s)S'^X{T/N). 
Corollary 2.4. Given the hypothesis of Lemma \2.3\ 

X{T/N)=SB-^S'^B where B^X{T/2N). 

Several more remarks about these factorizations are needed here. 

(1) In the case of time-reversing symmetry matrix, S is typically block diagonal with 
two blocks of opposite sign, one for the position variable and one for the momenta, 
that is, 

' F " 

-F 

where F is orthogonal. A matrix of this form is orthogonal and anti-commutes 
with J. 

(2) A matrix satisfying properties 3 and 4 of Lemma lZJl is symplectic with a multiplier 
of - 1 since S'^JS = -S'^SJ = J. 
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(3) If Y{s) is the fundamental matrix solution to (|2]i, then a similar argument shows 
that Y{-s + T/N) = SY{s)Yq^S'^Y{T/N) and consequently 

Y{T/N)=SYqB-^S'^B, where B = Y{T/2N). 

Applying this factorization theory results in expressing the matrix Yf^^Y{T), which is 
similar to X{T), as for some positive integer k, where the symplectic matrix W is the 
product of two involutions. If an eigenvalue of W lies on the unit circle, then so does its 

power The symplectic matrix W is called stable if all of its eigenvalues lie on the unit 
circle. 

Lemma 2.5. For a symplectic matrix W, suppose there is a matrix K such that 





(3) \{W + W-') = 



K 



Then W is stable if and only if all of the eigenvalues ofK are real and have absolute value 
smaller than or equal to 1. 

We will show for each of the periodic orbits under consideration, there is a choice of Yq 
such that W satisfies Lemma 12.51 This reduces the linear stability to the computation of 
the eigenvalues of K. As one of the eigenvalues of K is known to be real and have absolute 
value 1, the linear stability is determined by the numerical computation of one real number 
and showing that, within error, it lies between — 1 and 1 . 

3. Linear Stability for the Collinear Four-Body Symmetric Periodic 

Orbit 

The existence of the Schubart-like periodic orbit in the collinear four-body problem has 
been shown in Q. We review it here. For x\>X2> 0, we assume that four masses are 
located at x\, X2, —X2 and —xi with masses 1, m, m, and 1 respectively with m > 0. We 
also assume that the system remains symmetrically distributed about the center of mass 
located at the origin. The respective velocities of the four bodies are ii, X2,—X2,—xi where 
' = d/dt. The Netowanian equations are 

1 m m 

Ax\ {x\+X2f- [x\—X2Y 

m 1 1 

{^l+^lY {x\-X2Y 
We recount Sweatman's approach in lfT2l and lfT3l to regularize this system. The Hamil- 
tonian for this system is 

I 2 1 2 1 2»j 2m 

H^-;Wi + -7-W2 - :s ^ \ ' 

4 ^m lx\ 2x2 x\+X2 x\—X2 

where w\ = 2x\ and W2 = 2mx2 are the conjugate momenta to x\ and X2. Introduce new 
canonical coordinates qi,q2,p\,p2 by 

qi^xi~X2, q2 = 2x2, Pi=Wi, P2 = ^(wi+W2). 
The Hamiltonian in the new canonical coordinates is 

ff^}_fl^]_\p2 PiP2 , P2 2ot ot2 2m 1 



4\ mj m m q\ q2 qi+q2 '2.q\+q2 
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To regularize the equations of motion, Sweatman introduced a Levi-Civita type of canoni- 
cal transformation 

Qj^qi, Pi = 2Qipi (/=1,2), 

for the the canonical coordinates Qi , Q2,Pi,P2, and then replaced time t by the new inde- 
pendent variable s given by 

dt 



ds 



q\qI 



In the extended phase space, this produces the regularized Hamiltonian 



dt 
ds 



We fix the energy E 
(4) 



,2d2 



-Q\Q2P\P2 + Q\Pl 



Am 



1 . The Hamiltonian system in the new coordinate system is 



2 \ m m 



(5) 



Q'2 = ^ 

Zm 



Q1P2 - 2G2P1 



(6) P[ 



Am 



{Q2Pi-2QyP2)+2m^Q, 



AmQxQ% 



221 el 



2^2 



-2eie^, 



(7) 



Px 



Q1P2 Q2P1 1 



-4m22 + 



22^22, 



where ' is the derivative with respect to s. 

From the proof Q of the existence of the Schubart-like periodic orbit Qi{s), Q2{s), 
P\ {s), P2{s) of periodic T, there is a positive constant R{m) such that 

ei(0) =R{m), 22(0) =0, Pi(0) =0, P2(0) =2m3/2. 

These initial conditions correspond to a binary collision of the two inner bodies. By the 
construction of the periodic orbit, another binary collision of the two inner bodies occurs 
at i = r/2 where the conditions are 

ei(r/2) = -R{,n), e2(r/2) = 0, Pi(r/2) = 0, P2(r/2) = -2,n^^\ 

Simultaneous binary collisions correspond to the conditions of the periodic solution when 
s = T /4 and s ~ 3T /4, i.e., Qi (s) = at these values of s. The value of R{1) is approxi- 
mately 2.29559. Figure 1 contains a plot of the coordinates 2i , 22,^*1 ,^"2 of the periodic 
orbit when m= I. 
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Figure 1: The periodic solution in the coordinate system Qi,Q2,Pi,P2 
when m= I. 



3.1. Stability Reductions using Symmetry. The Schubart-like periodic solution y(,s') = 
(Gi {s),Q2{s),Pi {s),P2(s)) with period T in the collinear problem has two time-reversing 
symmetries. For 



the matrix 



1 








-1 


F 








-F 



is orthogonal and symmetric: S^^ = = S. It is also an involution, i.e., S = I. Since 
Sy{—s + r) is a solution of (|4]i through (I?), and since this solution shares the same initial 
conditions as 7(5) at 5 = by T-periodicity of y, uniqueness of solutions implies that the 
matrix S satisfies 

j{-s + T)=Sj{s) for alii. 

Thus 5 is a time-reversing symmetry of 7(5). With A' = 1, conditions (2), (3), and (4) in 
Lemma l23] are satisfied, and so by Corollarv 12.41 the monodromy matrix for y satisfies 

(8) X{T) = SX{T/2)-^S^X{T/2) = SX{T/2)-^SX{T/2). 

Consequently, from the above equation and = /, 

[SX{T)f= [X{T/2y'^SX{T/2)] [X{T /2y^SX{T/2)] =/. 

Since —Sy{—s + T/2) is a solution of (HI through (|7]i, and as —Sy{T/2) is the same as y(0), 
uniqueness of solutions implies that the matrix —S satisfies 

y(-5 + r/2) = -5y(i) for all s. 

Thus —S is another time-reversing symmetry of y(.s'). For N ~2, conditions (2), (3), and 
(4) of Lemma |2.3| are satisfied, and so Corollarv l2.4| implies that 



(9) 



X{T/2) = SX{T/4)-^SX{T/4) 
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For 

B=X{T/A), 

combining equations (|8]l and (|9]l gives 

X{T) = [SB-^SBf 
With A = SB-^SB and D = B'^SB, then 

X{T) =A^ = {SDf, 

where = I and = /. The two time-reversing symmetries S and —S of y are both 
involutions, and together they generate a symmetry group for y. 

3.2. A Good Basis. We have reduced the stabiUty analysis to the first quarter of the peri- 
odic oribt. Let Y{s) be the fundamental matrix solution to the linearized equations about 
Schubart-like periodic orbit y(i) with arbitrary initial conditions Yq. Let 

B = y(r/4). 

By the third remark following Corollarv 12.41 the matrix Yq'^Y{T), which is similar to the 
monodromy matrix = Y{T)Yq\ satisfies 

Yq-'Y{T) = {{Y^-'SYo) B-'SBf . 
The question of stability reduces to showing that the eigenvalues of 

W = {Yq^SYo)B-^SB 

are on the unit circle. An appropriate choice of Yq will simplify the factor Yq'^SYq in W. 
Set 



A = 



/ 

-/ 



Lemma 3.1. There exists Yq such that 

(1) Yq is orthogonal and symplectic, and 

(2) Yq-^SYq=A. 

Proof. Choose the third column of Yq to be y'(0)/||y'(0)|| = [0 1 0]^ = e2. For ea = 
[0 1 0]^,the matrix 

"0 10 

Yo^[Je2je3,e2,ej,]= o o Q 1 
-10 0^ 

is orthogonal and symplectic. Since S ~ diagjl, —1,-1,1}, it follows that Y^^SYq has the 
desired form. □ 

Setting D = B^SB and choosing Yq as constructed in Lemma lTTI gives 

w = {Yq ' syq) b-^sb = ad. 

The matrices A and D are both involutions, i.e., A^ = I, ~ I. From these it follows that 

W"' =DA. 

Because B is a symplectic matrix, a short computation using the formula for the inverse of 
a symplectic matrix shows that D has the form 

Li 

-L2 K 
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for 2x2 matrices K,L\ ,L2. It follows that 



and 



Hence, 



W 



I 

-/ 



K 



K 



K 



K 





1 01^. Setv 



Yq^Y{0). By the choice of Fq, 



We show that the first column of K is 

v^FoV(0) = ||y'(0)|k3. 
Since S is symmetric and Yq is orthogonal, then by the third remark after Corollary 12.41 



W = Y^;'SYqB-^SB = Yq'SYqB-'S'B^ Y^ Y{T/2). 

Now y'(i) is a solution of ^ = JD-r{y{s))t, and y'(0) = Y{Q)Yq^j'{Q) = Y{Q)v, and so 
y'(i) = Y{s)Yq^'y'{Q) = Y{s)v. This implies that 

Y^-W{T/2)^Y^Y{T/2)v^Wv. 

Since 7(5) satisfies y(-s + T/2) = -Sj{s) for all s, then y'{-s + T/2)= Sy'(s) for all s. 
Setting 5 = in this gives y'{T /2) = 5'y'(0). Since y'(0) is a nonzero scalar multiple of 62 
and since Se2 = —62, then 

Y,-WiT/2) = Y.-'siiO) = -Y.WiO) = -YoWiO) = v. 

Thus Wv = —V, implying that —1 is an eigenvalue of W and 63 is an eigenvector of W 
corresponding to this eigenvalue. Thus the first column of K is as claimed. The form of the 
rest of K comes from the formula for the inverse of a symplectic matrix and the definition 
ofD: 

-1 * 



K 







■ {SJC4) 



our numer- 



where c,- is the ;* column of Y{T /4-). 

3.3. Numerical Calculations. With an absolute error tolerance of 1 x 10"^', 
ical results for m= I showed that the initial condition 

ei(0)=/?(l)= 2.295592258717, 62(0) =0, Pi(0)=0, ^2(0)= 2 

leads to a periodic simultaneous binary collision periodic orbit (as in Figure 1) whose pe- 
riod T satisfies T/4 = 0.817348080989685. Using MATLAB and a Runge-Kutta-Fehlberg 
algorithm, we computed the columns of the matrix Y{T /4) with an absolute error tolerance 
of 4 X 10^^. From this, we got 

4 (57c4) =0.598490. 

For values of m between and 50 at 0.01 increments, we numerically computed the value 
of R{m) in the initial conditions and the value of the period T (with an absolute error 
tolerance of 4 x 10^^), and the values of C2 (57c4) (with an absolute error tolerance of 
1 X 10^^). The results of these computations are contained in Figure 2. 

A closer look at the numerical data in Figure 2 for where the value of 02(8 J 04) is close 
to 1 gives estimates of the two values of m where the stability of the periodic orbit changes. 
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The first critical value of m is approximately m = 2.83, and the second critical value of m 
is approximately m = 35.4. 

The eigenvalues of K are — I and Cj{SJc4). The eigenvalues of K are distinct for most 
values of m in [0,50] because of the rigorous numerical estimates we have for C2{SJc4.). 
Lemma |23] now implies the following linear stability result. 

Theorem 3.2. There exists small positive constants e,, / — 1,2,3,4 such that the periodic 
simultaneous binary collision orbit in the collinear symmetric four body problem with 
masses 1, m, m, 1 is linearly stable when m < 2.83 — Ei and 35.4 + 62 < m < 50, and is 
linearly unstable when 2.83 + £3 < m < 35.4 — £,4. 

This result confirms the linear stability analysis of Sweatman |[T3l for m between 
and 50, asserting that the periodic orbit is unstable when m is between 2.83 and 35.4. 
Simulations of the periodic orbit when m is between 2.83 and 35.4 indicate that the linear 
instability is manifested slowly over time. 




6 - 

0,5 - 

4 - 

3 - 



'0 5 10 15 2a 2S 3a 35 40 45 50 

Figure 2: The value of (5yc4) for values of m between and 50. 



4. Linear Stability for the 2D Symmetric Periodic Orbit 

In jS), we proved the existence of a special type of planar periodic solution of 2n bodies 
with equal masses. In this section, we are going to consider the linear stability of this 
periodic solution when n — 2. If (xi,X2) is the position of the first body, then the positions 
of the remaining three bodies are {x2,xi), {—xi,—X2), and (— X2,— xi). When each body 
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has mass m = 1, the Newtonian equations for this planar four-body problem are 

(jci -J[:2,X2 -xi) ixi,X2) (xi +X2,xi +X2) 



{■XUX2) 



2V2|xi-X2P 4(x2+x2)3/2 23/2|xi+X2P 



The initial conditions for the periodic orbit, and the periodic orbit are illustrated in Figure 
3. 




Figure 3: On the left are the initial conditions leading to the four-body 
two-dimensional periodic SBC oribt. On the right is the orbit. 

We adapt Sweatman's approach ( |[T2l . ifTSl ) to regularize this system. The Hamiltonian 
for this system is 

V2 ^/2 1 



H 



-W2) 



\Xl -X2\ 



\Xl +X2\ 



xj+xj 



where wi =4x1 and W2 = 4x2 are the conjugate momentum. In terms of the canonical 
coordinates {qi,q2,p\,p2) defined by 

^l=Xi-X2, ^2=Xi+X2, Wi=pi+p2, W2^P2-PI, 

the Hamiltonian becomes 



h^I{p\+pI) 



^/2 V2 



V2 



\qi\ 



1^2 1 



The Levi-Civita type of canonical transformation used to regularize the collinear problem 
now applies to the four body equal mass 2D problem. In terms of the canonical coordinates 
(G 1,62,^,^^2) defined by 

qi = Qf, Pi=^2QiPi (! = 1,2), 
and the new time variable s defined by 

the Hamiltonian in extended phase space becomes 

(10) F = ^(// -£) = UplQl+PlQ]) ^/2(e^ + G2) - -EQ]qI 



ds 
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where E is the total energy of the Hamiltonian H. The differential equations in terms of 
the new coordinates {2i,22,^'i,^'2} are 

(11) e'i = ^^ie2 



(12) 



(13) 



/'; = ->lGi+2x/2ei + ^^ 



^^^^^^ ■ 



1 9 r 2V2Q2Q] 

(14) P'2 = --PlQ2+2^/2Q2 + ^^^^ 



"^2^1 '2EQ2Q\. 



Unlike the collinear problem, we do not fix the value of E here. As shown in Q, for each 
^ > there exists vo > such that the initial conditions 



(15) 



gi (0) = C, G2(0) = C, ^i (0) = -4vo, ^^2(0) = 4vo, 



lead to a periodic solution with a minimal period T . From F = 0, the value of E is deter- 
mined by this choice of ^ and vq. By its construction in fSl, this periodic orbit satisfies 



ei(r/4) = -C, e2(r/4) = C, Pi(r/4) = -4vo, P2(r/4) 



-4vo. 



Simultaneous binary collisions correspond to i = r/8,5r/8 i.e., when 2i('5) = 0. and to 
5 = 3r/8,7r/8, i.e., when 22(^) = 0. For ^ = 1, 4vo = 2.57486992651942, and r/8 = 
1.62047369909693. Figure 4 illustrates the coordinates (Gi,G2,^'i,^'2) of this periodic 
solution. 




Figure 4: The periodic solution in the coordinate system Gi , G2,^'i ,^'2- 
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4. 1 . Stability Reductions using Symmetry. We will reduce the stability analysis to the 
first eighth of the periodic orbit. The symmetric periodic 2D orbit 

i{t)^{Qi{t)M)Mt)Mt)) 

with period T has a time-reversing symmetry and a time-preserving symmetry. For 



the matrices 



Sf = 



-1 

1 



F 
F 



G = 



Sg = 



-1 




satisfy Sp^ = 5^, ^ I, Sj, ^ /, Sj- = /, 4 = = ^G, and (5f5g)^ = /. Since 
Y(i + r/4) and Sfjis) = i-Q2{s),Qi{s),-P2is),Pi{s)) are solutions of ([TTTi through ([T4li 
and share the same initial conditions when s = 0, uniqueness of solutions implies that 



Y(i + r/4) = 5f y(i) for all s. 

Thus Sf is a time-preserving symmetry of y(s). With N = 4, conditions (2), (3), and (4) of 
Lemma IZTI are satisfied, so that Corollarv l2.2l (with A; = 4) and Sp^I imply that 

X{T) = St {slx{T/4-)Y = {S^pX{T/4)Y . 

Since y(— i + r/4) and Sajis) are solutions of (fTTl i through (fl4] i and share the same initial 
conditions when s = 0, uniqueness of solutions implies that 

yi-s + r/4) = SgJ{s) for all s. 

Thus Sg is a time-reversing symmetry for 7(5). With N = 4, conditions (2), (3), and (4) of 
Lemma |23] are satisfied, and so Corollarv 12 .41 implies that 



X{T/4) ^ Sg [X{T/^)]-' 5jX(r/8) = Sg [X{T/S)]-' SgX{T/^). 



Let 



B=X{T/8). 



Combining the factorization ofX{T) that involves Sp and the factorization ofX{T/4) that 
involves Sg gives the factorization 



X{T)^{S''pSgB-^SgB)\ 



Setting 



Q = S'pSc 



1 

1 














-1 





and D = B ^SgB results in the factorization 



X{T) = (QD)'^ 

where Q and D are both involutions. The symmetries Sf and Sg generate a D4 symmetry 
group for the periodic orbit y(i). 
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4.2. A Good Basis. Let Y{s) be the fundamental matrix solution to the linearized equa- 
tions about the 2D periodic orbit ^{s) with arbitrary initial conditions Yq. Let 

B = y(r/8). 

By remarks following Corollaries 12.21 and l2 .41 the matrix Yq^Y(T), which is similar to the 
monodromy matrix = Y{T)Yq'^, satisfies 

Y^-'YiT) = (Yq-'sISgYoB-'SgB)'^ = {Yo-^QYoB-'SgB)\ 
The question of linear stability reduces to showing that the eigenvalues of 



W = Yq^QYoB-^SgB 



are on the unit circle. Recall that 



A = 




/ 



Lemma 4.1. There exists Yq such that 

(1) Yq is orthogonal and symplectic, and 



(2) 7-^270 = A- 
Proof. Choose the third column of Yq to be 



Y'(0) 



1 



— —a a 
c '■ 



llY'(0)|| 

where a = voC^/2, b = E^"' = (2vf) - 2^2 - 1 )C and c 
the /* column of Yd. Define 



b b\ 



\/2^ 



'2b^. Let col,(Fo) denote 



co\i{Y(i)=J ■coUYq) = -\ b b a -a] . 

We now choose col4(Fo) such that col4(}o) is orthogonal to col3(}o), and col4(}o) is 
one of the eigenvectors of Q with respect to its eigenvalue of — L Since the eigenspace of 
Q corresponding to its eigenvalue of — 1 is 

span{ [1 -1 [ 1 1 ]^}, 



define 



and 



The matrix 



coi4(yo) 



- \ b -b 

c 



col2(Fo) = ^-0014(70) = -[ fl a -b b]' 



Yo 





b 


a 


—a 


b 


1 


b 


a 


a 


-b 


c 


a 


-b 


b 


a 




—a 


b 


b 


a 



is both symplectic and orthogonal and it satisfies Yq ' QYq = A 



□ 



Setting D ~ B^^SgB and choosing Yq to be the matrix constructed in Lemma |4TT] gives 
W = AD. The matrices A and D are involutions (the latter because Sq ~ I). As in Section 
3.2, W^' = DA, and there is a 2 x 2 matrix K such that 



1 
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We show that the first column of K is [1 0]^. Since 
the third remark following Corollary 12 .41 that 

W = Yq^sISgYqB-^SgB = Yq^sIy{t/a) 

Set V = YQ^y'{Q). By the choice of the matrix Yq, 

v = J'o"V(0)=FoV(0)- 

Because y'{s) is a solution to the linearized equation 4 = JD~T{y{s) )i, and because y'(0) = 
Y{0)YQ^i{Q), then Y'(i) = Y{s)YQ^i{0) for all s. Hence, 

(16) Wv = Y^SlY{T/4)v = Y^Sly'{T/4). 

Since y satisfies y{s + T/4) = Spyis) for all s and 5^ ' = 5^, it then follows that 

y'{s)=Sp'y'{s + T/4)=sW{^ + T/4). 
Setting i = in this gives y'(0) — Spy'{T /4), and consequently that 

(17) Yjs'j:y'{T/4) = Yly'(O) = yo"'Y'(0) = v. 

Equations (fT6l ) and il% now combine to show that Wv — v, i.e, that 1 is an eigenvalue of 
W and 63 is an eigenvector for W corresponding to this eigenvalue. The first column of K 
is as claimed. The rest of K comes from the formula for the inverse of a symplectic matrix 
and the definition of D: 

r 1 * 

^ = [ 4iSGJc4) \ ' 

where c,- is the /* column of B = Y{T /Bi). 

4.3. Numerical Calculations. Having not fixed E, we used an invariant scaling of the 
coordinates and time in equations (fTTT i through (fl4] i to preselect a period T before numer- 
ically computing the initial conditions for a periodic simultaneous binary collision orbit. 
For e > 0, if (2i i^^)^ Qii^), Pi (s), Pii^) is a periodic simultaneous collision orbit of equa- 
tions ( fTTT i through (fl4] i. then replacing E with e^^E shows that egi (e^), eQ2{£s), Pi (es), 
P2{es) is also a periodic simultaneous binary collision orbit with energy e^'£ and period 
g-i J Furthermore, it is straight-forward to show that monodromy matrices for the peri- 
odic simultaneous binary collision orbits corresponding to values of e 7^ 1 are all similar to 
that for e = 1 . Thus the linear stability of a periodic simultaneous binary collision orbit for 
one e > implies the linear stability of the periodic simultaneous binary collision orbits 
for all e > 0. 

We rigorously computed the value of Cj{SgJc4) for the periodic simultaneous binary 
collision orbit whose period is T = 8. This means that the first time of a simultaneous bi- 
nary colUsion for this orbit is at = 1 . We set Qi{0) = Q2 (0) = ^ and -Pi (0) = ^2(0) = r|, 
and defined a function F(^,r|) to be equal to the vector quantity (Qi (1)5^2(1))- We used 
Newton's method and a good initial guess to find a root (^,r|) of F. This involved com- 
puting the Jacobian of F which was done using the linearized equations. With an absolute 
error tolerance of 6 x 10^ this numerical method shows that the initial conditions 



= Sg, = ^0^' it follows by 
= yJs]:y{T/4). 




-||Y'(0)lk3. 



Gi(0) = 62(0) = 1.62047369909693, -Pi (0) = ^2(0) = 2.57486992651942, 
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lead to a periodic solution with a period of T = 8, and a value of £ « —1.142329388. 
Using MATLAB and a Runge-Kutta-Fehlberg algorithm, we computed the columns of the 
matrix 7(7/8) with an absolute error tolerance of 2.5 x 10^'^. From this we got 

4 (SgJca) = -0.68024151010592. 

Using the scaling of coordinates and time described above, the initial conditions for the 
periodic simultaneous binary collision orbit shown in Figures 3 and 4 are 

Gi(0) =62(0) = 1, -Pi (0) =P2(0) =2.57486992651942 

with a period T satisfying r/8 = 1.62047369909693, and energy E « -2.999682732. 

For the periodic simultaneous binary collision orbit, the rigorous estimate of the eigen- 
value C2 (5'g/c4) of K and its distinctiveness from the eigenvalue 1 of AT combine with 
Lemma |23] to give the following stability result. 

Theorem 4.2. The periodic simultaneous binary collision orbit in the ID-symmetric equal 
mass four-body problem is linearly stable. 

When c^(5g/c4) is real and between —1 and 1, it is the real part of an eigenvalue 
with unit modulus for W (see Q)- For the periodic simultaneous collision orbit, the real 
part of exp(37t//4), that is —(1 /2)a/2, is fairly close to the rigorously estimated value of 

{SgJca)- Raising exp(37t//4) to the fourth power gives exp(37t;) = — 1, and so two of 
the eigenvalues of the monodromy matrix of the periodic simultaneous binary collision 
orbit are close to —1. The symmetry reductions used to compute the eigenvalues over just 
one-eighth of the period and the rigorous estimate of c^ {SgJca) showing that it is clearly 
between — 1 and 1, assures the linear stability of the periodic simultaneous binary collision 
orbit. 
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